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An analytical method to calculate the rotor excitation forces arising from 
labyrinth seals is presented. The objective is to model the gas flow through the 
seal clearance passages and cavities when the rotor is positioned eccentric relative 
to the stator center. The seal flow model used in the analysis yields solutions 
which validate the experimentally observed influence of seal parameters on seal 
forces reported in the literature. The analytically predicted seal pressure 
distributions and forces have been compared with published experimental results. 


INTRODUCTION 

Whi le performing the function of sealing, labyrinth seals have been identified 
as a potential source of destabilizing forces which cause self excited rotor whirl 
of turborotors. Typically these seals are employed to effect sealing of the working 
fluid between two stages of a compressor or turbine and in adjacent cooling air 
compartments which are maintained at different pressures. Labyrinth seals which 
have a configuration as shown in figure 1 are designed to minimize the leakage flow 
by successively expanding the gas past the or if ice- 1 ike clearance passages at the 
seal fin locations. 

An analytical model describing the excitation of turborotors by labyrinth seals 
was first published by Alford (ref. 1). In the recent literature (ref. 2 to 7) 
experimental and analytical methods of predicting labyrinth seal forces have been 
reported. Experimental results presented by Benckert and Yfachter (ref. 2) show that 
a restoring force along the line of centers and a lateral force normal to the line 
act on the rotor when the rotor center is displaced to an eccentric position with 
respect to the seal center as shown in figure 2. The present work is aimed at 
analytically predicting these seal forces which result from the circumferentially 
varying pressure distributions within the seal cavities. 


SYMBOLS 


c 

seal radial clearance 

Cc 

orifice flow coefficent 

Ce 

energy recovery factor 

Co 

nominal radial clearance 

e 

rotor eccentricity 

®r 

eccentricity ratio (e/Co) 

f 

circumferential flow area 

Fr, Ft 

radial and tangential seal 

h 

seal fin height 
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Subscripts 

i 

j 

r 

s 


cavi ty width 

number of circumferential nodes 
rotor rpm 

number of seal fins 

pressure (Po- seal inlet, Pe- seal exit) 

axial mass flow rate per unit circumferential length 

gas constant 

reynolds number 

rotor radius 

axial pitch 

gas temperature 

circumferential flow velocity 

inlet swirl velocity 

length of wetted surface 

flow coefficient 

friction coefficient at the walls 
coefficient of viscosity 
reference angle shown in figure 2 
gas density 
shear stress 

circumferential coordinate 
angular velocity of rotor rotation 

fin and cavity number designations 
circumferential node number 
rotor 
stator 


GOVERNING EQUATIONS 


asymmetry caused 
fluid mechanics 


The intricate flow path of labyrinth seals, combined with the 
by rotor eccentricity complicate the analysis employing fundamental 
governing equations. A relatively simpler set of labyrinth seal governing equations 
derived by Kostyuk (ref. 8) are employed in the analysis here. For the fluid 
macro-element shown in figure 1, the following conservation of 
circumferential momentum equations are derived, neglecting the 
small clearance-to-radius ratios. 


mass and 
curvature 


effect for 


2 (pjfj) + 1 


3t 


fts 30 


(pjUjf j ) + (q i+1 - qj ) =0 


( 1 ! 


2 (pjUjfj) + 2 _ L (Pj u 2 | f | ) + (q i+1 Uj - q j u j _i ) 
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( 2 ) 


Frictional resistance at the seal walls are calculated using the following 
re lat i onshi ps. 
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(3) 


si 


= *si W 


r i 


= Vi 


p: (wRs-u; ) 2 


where: 


A = _£ 0<Re<1000, A = °- 045 1000 < Re 

Re Re- 25 

The Reynolds numbers are given by: 


u i| h iPi Rp ; = |^S-Uj 
~ M M 2p 

Proper signs are attached to the friction coefficient A. based on the relative 
velocity of the fluid with respect to the rotor and stator walls. Assuming 
isothermal flow process, Martin's (ref. 9) leakage flow equation for small pressure 
drop accross each fin is given by: 

qj = at j (P 2 j_i - P 2 j ) ' ® (5) 



|Mi 


(4) 


where: 


a: = 


1 Ccj 

/RT ( 1-Ce j ) 6 


( 6 ) 


Orifice flow coefficients Cc for annular orifices presented by Vermes (ref. 10) 
have been used in the calculations. The residual energy factor Ce accounts for the 
part of the kinetic energy in the leakage flow from the (i-l)th orifice that is 
present at the inlet to the i-th orifice. Representing the orifice flow by one half 
of a flat symmetrical jet. Vermes (ref. 10) derived the following expression for the 
residual energy factor: 


Ce- = u.-Ji. 

1 (I j/Cj) + 7.23 (7 ) 

Equations (1), (2), (5) and the ideal gas relationship 

Pj = pj RT (8) 

form the system of equations to be solved for the circumferential pressure 
distribution resulting due to the rotor eccentricity. Although figure 1 shows fins 
on the rotor, the equations derived are also val.id for a labyrinth with fins on the 
stator . 


NUMERICAL PROCEDURE 

For the rotor undergoing steady rotational motion about a given eccentric 
position (fig. 2), we can drop off the time derivative terms in equations (1) and 
(2). After substituting for qj and p j from equations (5) and (8), equations (1) and 
( 2 ) become : 
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( 10 ) 


The seal clearance and circumferential flow area for the eccentric rotor 
configuration are given by: 

Cj = Co; - e cos0 (11) 

fj = lj (hj + Coj - e cos0) (12) 

Equations (9) and (10) are nonlinear ordinary differential equations to be 
solved to get the pressure and velocity distributions, Pj (0) and uj (0) in each 
seal cavity. 


After introducing the central finite difference approximations given below we 
divide the circumferential coordinate length, Rs0 into m equal parts. 
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_ p i j+1 ~ p i j-1 
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The resulting system of nonlinear algebraic equations are represented 
symbo I i ca I I y as : 

' F (Pj j , uj j ) = 0 (15) 

G ( p j j, Uj j) = 0 i=1,2..(n-1) j=1,2..m (16) 

Equations (15) and (16) are solved employing the Newton Raphson iteration 
procedure in the multidimensional space formed by the variables Pjj and Ujj. After 
dropping the second and higher order terms from the Taylor series expansions of 
equations (15) and (16), we arrive at the iteration scheme given by: 
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By assuming an initial guess for the vector [Pjj Ujj]^, the Jacobian matrix and 
the function values are evaluated to calculate the new vector of variables 
[Pjj Ujj]jj + i- The iterations are continued until the change in each component of 
the solution vector is less than a predetermined error value. The seal forces are 
obtained by integration of the pressure distribution over the rotor surface using 
the summation equations given by: 


n i m 

Fr = RsA0 j. lj £ Pjj cos0j 

i-1 j-1‘ 


Ft = RsAd I j £ Pjj sinOj 
i =1 j=1 


(18) 

(19) 


RESISTS 

Seal cavity pressure distributions obtained were compared with the experimental 
results of Leong and Brown (ref. 4) shown in figure 3. The 11-cavity seal having 
fins on the stator with inlet and exit pressure of Po=2.914 bar, Pe=1 bar was 
investigated at a rotor speed of 7504 rpm and eccentricity ratio of 0.6. The 
analytically predicted pressure distribution patterns in the first two seal cavities 
are in agreement with the experimentally measured patterns, with a peak pressure 
occuring near the 180 degrees location in the first cavity while the pressure is 
minimum near this location in the second cavity. Pressure levels predicted by the 
theory beyond the first three cavities are in fair agreement with the experimental 
results. The difference in cavity pressure levels between the theory and experiment 
is attributable to possible difference between the orifice flow coefficients used in 
the calculations and the actual flow coefficients in the experiment. Also, possible 
random variations of the seal clearance at successive fin locations in the test seal 
could contribute to the difference between the experimental and analytical results. 

The experimental results of Benckert and Wachter (ref. 2) were used to verify 
the seal forces calculated analytically. Figures 5 and 6 show the variation of seal 
forces with eccentricity ratio for a 17-cavity seal with fins on the stator. 
Tangential forces generated when the rotor rotates at 9482 rpm and 4773 rpm, for the 
seal inlet and exit pressures of Po=1 bar, Pes" 1 bar are shown in figure 4(a). The 
negative eccentricity ratios indicate that the rotor center has been displaced along 
the negative x-axis in figure 2. The analytically predicted tangential seal forces 
are in good agreement with the experimental results showing an average deviation of 
about 16 percent, excluding the zero eccentricity ratio case where the experimental 
results show a nonzero force. Radial forces predicted by the analysis are shown in 
figure 4(b); experimental results on radial forces were not reported in the 
reference for the case investigated. These results show that the seal force 
magnitudes increase with rotor rotational speed for any given eccentricity ratio. 

The tangential forces act in a direction that would initiate or promote forward 
whirl motion of the rotor center. The radial forces are decentering where they tend 
to push the rotor away from the seal center. 

The seal forces variation with eccentricity ratio for zero rotor speed, but 
nonzero swirl velocity at the seal inlet are presented in figures 5(a) and 5(b). 

The seal pressures in this case are Po=2.041 bar, Pe=1 bar. The tangential seal 
forces have been predicted quite close to the experimental results, the maximum 
difference being 17 percent. The radial forces from the experiment, again were not 
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reported in the reference. The two different inlet swirl velocities investigated 
snow that the seal forces increase with the inlet swirl velocity and appreciable 
magnitudes of tangential forces are generated which would promote forward whirl 
instability in rotor systems. The radial forces in this case are of relatively 
small magnitude and they act in a centering direction tending to restore the rotor 
center back to the seal center. 


CONCLUSIONS 

Employing the seal flow governing equations derived by Kostyuk (ref. 8), 
methods of predicting seal forces for an infinitesimal eccentricity have been 
reported in the literature where seal stiffness coefficients are presented. In this 
analysis the radial and tangential components of seal forces for large eccentricity 
ratios have been predicted, for the first time, employing Kostyuk's equations. The 
analysis yields results which validate the experimentally observed variation of seal 
forces with rotor eccentricity, inlet swirl velocity and rotor rotational speed. 

The tangential force components predicted by the theory, especially the swirl 
flow-induced forces are in good quantitative agreement with the experimental 
results. The analysis also shows the nonlinear dependence of the forces with 
eccentricity beyond e^ .4. This behavior needs to be given appropriate 
consideration in large amplitude rotor dynamic investigations where nonlinear 
stability problems may be encountered. 
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Labyrinth seal fluid element 
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Figure 3. - Seal pressure distribution. 
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Figure 4. - Seal forces with variable rotor speed. 
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(a) Tangential force Ft. 

(b) Radial force Fr. 

Figure 5, - Seal forces with variable inlet swirl velocity. 
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